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Fermi Large Area Telescope (LAT) has recently detected 8 gamma-ray millisecond pulsars 
(MSPs), providing an unprecedented opportunity to probe the magnetospheres of these low- 
spin-down pulsars. We performed 3D emission modeling, including various Special Relativistic 
effects, in the context of pair-starved polar cap (PSPC), slot gap (SG), and outer gap (OG) 
pulsar models. Most of the light curves are best fit by SG and OG models, surprisingly 
indicating the presence of narrow accelerating gaps limited by robust pair production. All 
model fits imply high-altitude emission, and we observe exclusive differentiation of the current 
gamma-ray MSP population into two sub-classes: light curve shapes and lags across wavebands 
impose either PSPC or SG / OG-type geometries. 



1 Introduction 

Millisecond pulsars (MSPs) are rapidly-rotating neutron stars, characterized by relatively short 
spin periods and low surface magnetic fields. It is thought that transfer of mass and angular 
momentum from a binary companion during an accretion phase, may revive their "dead" (radio- 
invisible) younger pulsar progenitors which have "spun down" previously. ^ MSPs are thus 
believed to be "recycled" pulsars,^ decending from lo w-mass X-ray binaries (LMXRBs).^ 

MSPs were predicted to be visible in gamma ra y S- l3ISIHIZIHMl0] Observational confirmation 
was provided by the Fermi Large Area Telescope (LAT), a high-energy (HE) gamma-ray tele- 
scope with a large field of view (2.4 sr), operating in the energy range ~ 20 MeV to > 300 GeVUJJ, 



which recently detected 8 gamma-ray MSPs. 1 ^ 

Mainly two different models have been employed to describe HE radiation from pulsars. 
Polar cap (PC) models^] assume injection of primary electrons from the stellar surface. These 
primaries radiate HE curvature radiation (CR) or inverse Compton scattering (ICS) gamma rays 
which are converted into electron-positron pairs via magnetic pair production. A low-altitude 
pair formation front (PFF) is established close to the surface, screening the accelerating electric 
field. El By allowing variation of the CR PFF altitude across the PC, a slot gap (SG) is formed 
along the last open magnetic field lines. This gap extends from the neutron star surf ace to near 
the so-called light cylinder, and allows acceleration of primaries up to high altitudes. DSHD The 
SG mo del is a p ossible physical realization of the two-pole caustic geometry.^3 Outer gap (OG) 
models D-8 | 19 | 20 | assume that HE radiation is generated by cascades of electron-positron pairs 
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Figure 1: Left: Phaseplots for PSR J0030+0451. Panel (a) is for an SG model with (a,Q = (70°, 80°), while 
panel (b) is for an OG model with (a, <") = (80°, 70°). Right: Observed and fitted light curves for PSR J0030+0451. 
In panel (a), we show the Fermi-LAT data (histogram) with estimated background level (dashed line), and SG 
(magenta line) and OG (green line) fits. Panel (b) shows the radio data (blue line), and the magenta and green 
lines correspond to the same (a, C) combinations as those of the SG and OG fits in panel (a). 



produced via photon-photon pair creation. This emission occurs close to the last open field lines, 
above the "null-charge surface", and is thus confined to the outer magnetosphere. Recently, a 
3D OG solution was found^^ which extends toward the NS surface. For M SPs with mostly 
unscreened magnetospheres, a "pair-starved polar cap" (PSPC) mo deli™ applies in which 
charges are accelerated up to high altitudes over the full open-field-line region. 

We study the Fermi-LAT MSP population^^ using 3D emission modeling, including Special 
Relativist ic (SR) e ffects of aberration and time-of-flight delays, and rotational sweepback of B- 
field lines ESEIEE] [ n the context of geometric SG and OG models, as well as a CR PSPC model, 
and obtain fits for gamma-ray and radio light curves. Our calculations are complementary to 
another study^that focuses on younger pulsars. More details are provided elsewhere. 

2 Model Description 

We use an implementation ^® ! ^ of the retarded vacuum dipolar B-field solutionis G f a rotator 
inclined by an angle a with respect to the rotation axis. The PC shape is distorted asymmet- 
rically by rotational sweepback of field lines. We calculate the curvature radius of the B-field 
lines in the inertial observer frame, and assume SG and OG gaps to be confined between two 
field lines with footpoints close to the PC rim, and having constant emissivity over the emit- 
ting volume. In the PSPC case, we collect photons with energies above 100 MeV from the full 
open-field-line region. 

For the PSPC regime, we only c onsider CR losses suffered by electron pri maries moving 
along the field lines. Previous 

studies lMZEEQ] 

have used lower-altitude 

solutionsEEQ] for E 

the PSPC E-field parallel to the B-field. As Fermi-LAT results seem to indicate that the HE 
radiation is originating in the outer magnetosphere,^] we now incorporate the solutionis for 
altitudes close to the light cylinder in the small-angle approximation by matching this solution to 
the two lower-altitude ones. We find conservation of energy when solving the transport equation 
for relativistic electron primaries using the full E-field solution. 

We use an empirical radio cone model based on the characterization^^ Q f the radio emission 
as magnetic axis-centered core and conal beams. We adopt a description^ZI based on fits^S of 
average-pulse profiles of a small collection of pulsars at 400 MHz. The conal emission occurs 



Table 1: Model fits for inclination and observer angles a, £. 
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at altitudes of 10% — 20% of the light cylinder radius. We lastly use relative units for the cone 
beam luminosity. 

3 Results 

SR effects of aberration and time-of-flight delays cause phase shifts that nearly cancel those 
due to the curvature of the B-field on trailing field lines, leading to accumulation of emission 
(initially emitted tangentially to the B-field lines) around narrow phase bands 4>. This results 
in caustic structures (on phaseplots of C vs. 0, with C the observer angle) around <fi ~ 0.0 — 0.1 
and 4> ~ 0.4 — 0.6 in phase for SG and OG models. No emission originates below the null charge 
surface in the OG model, so that an observer can only see emission from one magnetic pole, in 
contrast to the SG models. 

We studied a large number of light curves for each of the different pulsar models. Both OG 
and SG models have a preponderance of double-peaked light curves at similar phases. OG models 
do not exist at all angle combinations, while SG models do (due to emission occuring below the 
null charge surface). One may find sharp, solitary peaks for some regions in (a, C)-space in OG 
models, while the corresponding SG-peaks usually have additional low-level features. The PSPC 
model have mostly single-peaked gamma-ray profiles roughly in phase with the (single) radio 
peaks, but one may see radiation from both poles for large a, leading to double-peaked profiles 
separated by ~ 0.5 in phase. The radio peak multiplicity depends on the observer's geometry. In 
significantly off-beam geometries, only the relatively larger gamma-ray cone is seen, explaining 
the phenomenon of "radio-quiet" pulsars. 

We chose best-fit light curves from the various models to match the MSP gamma-ray and 
radio data by eye. As an example, the left panel of Figure Q] shows phaseplots associated with 
the best light curve fits (right panel) obtained for PSR J0030+0451. We inferred values for a 
and C for each MSP (Tabled]), and compared them with (somewhat uncertain) values obtained 
from radio polarimetric measurements (references for radio data elsewhere^l) . 

4 Discussion and Conclusions 

We have compared 3D model predictions of gamma-ray and radio radiation with MSP gamma- 
ray data from Fermi-LAT in the framework of geometric SG and OG pulsar models, and also for 
the full-radiation PSPC model. Surprisingly, some MSPs have double-peaked light curves well 
fit by SG / OG models, indicating strong screening of Eu by pair production. The larger radio 
beam widths of MSPs compared to those of canonical pulsars furthermore implies relatively few 



radio-quiet MSPs. We found exclusive differentiation between the SG / OG models on the one 
hand, and the PSPC model on the other hand. Emission in all models considered comes from 
the outer magnetosphere. Our fits of a and ( are in reasonable agreement with values inferred 
from MSP radio polarization measurements. For PSR J0437-4715 and PSR J0613-0200, the SG 
model predicts a small precursor to the main gamma-ray peak, but not the OG model. This 
may become a future model discriminator. Future phase-resolved spectroscopy made possible 
by the quality of Fermi-LAT data should challenge pulsar models. 
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